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^ : ABSTRACT 
> ■ 

o : 

^ I We present generalized supernova (SN) light curve (LC) models for a variety of 

^ I power inputs including the previously proposed ideas of radioactive decay of ^^Ni 

I and ^^Co and magnetar spin-down. We extend those solutions to include finite 

pq I progenitor radius and stationary photospheres as might be the case for SN that 

[nr] I are powered by interaction of the ejecta with circumstellar matter (CSM). We 

I provide an expression for the power input that is produced by self-similar forward 

^1 and reverse shocks that efficiently convert their kinetic energy into radiation. We 

O I find that this ejecta-CSM interaction luminosity that we derive is in agreement 

^ \ with results from multi-dimensional radiation hydrodynamics simulations in the 

1 case of an optically-thin CSM. We develop a semi-analytical model for the case of 
^ I an optically-thick CSM by invoking an approximation for the effects of radiative 
I diffusion similar to that adopted by Arnett (1982) for SN II and compare this 

^ I model to the results of numerical radiation hydrodynamics models. This model 

^ I can give complex light curves, but for monotonically declining shock input, the 

LCs have a smooth rise, peak and decline. In the context of this model, we 
I I provide predictions of the shock breakout of the forward shock from the optically- 

^ ' thick part of the CSM envelope. We also introduce a hybrid LC model that 

I incorporates ejecta-CSM interaction plus ^^Ni and ^^Co radioactive decay input. 

' We fit this hybrid model to the LC of the Super-Luminous Supernova (SLSN) 

' 2006gy. We find that shock heating produced by ejecta-CSM interaction plus 

■ - - ' some contribution from radioactive decay provides a better fit to the LC of this 

event than previously presented models. We also address the relation between 
SN IIL and SN Iln with ejecta-CSM interaction models. The faster decline of SN 
IIL can be reproduced by the diffusion of previously deposited shock power if the 
shock power input to the diffusive component vanishes when the reverse shock 
sweeps up the whole ejecta and/or the forward shock propagates through the 
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optically-thick CS matter. A CSM interaction with forward and reverse shock 
power input can produce the LCs of SN Iln in terms of duration, shape and 
decline rate, depending on the properties of the CSM envelope and the progenitor 
star. This model can also produce LCs that arc symmetric in shape around peak 
luminosity, which is the case for the observed LCs of some recently discovered 
peculiar transient events. We conclude that the observed LC variety of SN Iln 
and of some SLSNe is hkely to be a byproduct of the large range of conditions 
relevant to significant ejecta-CSM interaction as a power source. 

Subject headings: circumstellar matter - stars: light curve - supernovae: general 
- supernovae: individual - light curve 



1. INTRODUCTION 

The current SN classification scheme is based on both the properties of the spectrum 
and of the light curve. The basic properties of stripped envelope SNe (Ia,b,c) hght curves 
(LCs) are well reproduced by considering the diffusion of the radioactive decay energy of 
^^Ni and ^^Co into homologously expanding SN ejecta (Arnett 1979, 1980, 1982, 1996 - 
hereafter A79, A80, A82, A96). Core-collapse SNe are divided into the Type Ib/c and Type 
II subclasses. SN Ib/Ic are considered explosions of compact progenitors that have lost their 
outer hydrogen (for SN lb) and helium (for SN Ic) envelopes. A79 addressed what we would 
now call SN lb; A80 and A82 pertained more directly to SN la, but the principles are the 
same. For Type II SN, the following subtypes have been proposed: Type IIP, Type lib. 
Type IIL and Type Iln. SN IIP explosions are the most common, and they are beheved 
to be the result of the death of a massive red supergiant (RSG) progenitor star. The LCs 
of SN IIP are characterized by a long plateau that is indicative of the energy hberated by 
a recombination front in the extended hydrogen envelope of the progenitor star. The late 
time decline rate of many SN IIP LCs is consistent with that of the radiactive decay of 
^^Co. Arnett & Fu (1989) presented a semi- analytical model that incorporates the effects 
of radiative diffusion with energy deposition from the radioactive decay of ^^Co (also pulsar 
and fallback energy) plus H recombination in order to fit the observed LC of SN 1987A. This 
model can produce a variety of SN IlP-like LCs (and probably SN lib, although this has not 
been explored in depth) depending on the choice of H mass and the composition of the outer 
shell of the progenitor RSG star. More accurate radiation hydrodynamics simulations of SN 
IIP have been done yielding simultaneous LCs and spectra (sec examples in Falk & Arnett 
1977; Falk 1978; Klein & Chevalier 1978; Kasen & Woosley 2009; Dessart, Livne & Waldman 
2010; Bersten et al. 2011; Dessart & Hiller 2011) as well as models of SN 1987A that have 
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reproduced the observables in great detail (Blinnikov et al. 2000; Dessart & Killer 2010). 
One- dimensional numerical LC models based on radiation hydrodynamics calculations are 
available for SN IIL (Swartz ct al. 1991; Blinnikov & Bartunov 1993), but a simple semi- 
analytical approach, similar to that presented by A82, A96 for the case of SN la has not 
been provided. The simulations by Swartz ct al. (1991) were able to reproduce some of the 
characteristics of SN IIL LCs under the assumption of 0-Ne-Mg (Oxygen-Neon-Magnesium) 
core-collapse, but the fits were not good for all phases of the observed LCs. The LCs of SN 
Iln present a diversity when it comes to shape, peak luminosity, duration and decline rate. 
Some of the recently discovered Super-Luminous Supernovae (SLSNe) are classified as SN 
Iln events, mainly due to the spectroscopic signatures of ejecta-CSM interaction (Schlegel 
1996). 

Although there are numerical simulations of the interaction of SN ejecta with the CSM 
and model LCs under optically-thin conditions (Falk & Arnett 1977; Van Marie et al. 2010) 
the only attempt made so far to incorporate the effects of radiative diffusion of shock- 
deposited energy in an optically-thick CSM envelope was done in the case of red supergiant 
(RSG) progenitors surrounded by a CSM with a limited set of density profiles {pcsM ~ 
j,-2,-i.5^ (Moriya et al. 2011). The work done by Moriya et al. (2011) using the 1-D 
radiation hydrodynamics code STELLA (Blinnikov & Bartunov 1993; Blinnikov et al. 1998; 
Bhnnikov et al. 2006) reproduced the LC of ultraviolet-rich SN IIP SN 2009kf reasonably well 
and provided a tool to study the LCs of low-luminosity SN Iln powered by the interaction of 
the SN ejecta with a CS wind. For higher luminosity events in the regime of SLSNe, Moriya 
et al. (2011) find that extraordinary mass-loss rates are required that are inconsistent with 
an RSG progenitor. This and other factors motivate us to update the analytical models of 
A80, 82 and determine approximate LCs for SLSNe powered by optically-thick SN ejecta 
- CSM interaction for a variety of CSM and progenitor characteristics. This is important 
in order to produce a physically acceptable model for a SLSN Iln LC for which there is 
a smooth rise to maximum light; a simplified model of instantaneous shock heating yields 
a LC that increases steeply to maximum in a very short time-scale. The large parameter 
space that is connected with SN Iln (properties of the SN progenitor, the SN ejecta and of 
the CSM envelope) provide a natural explanation for the observed diversity of SN Iln LCs. 
Another motivation for developing an efficient approximate analytical model is to perform a 
qualitative study of this large parameter space. Given that many of the recently discovered 
SLSNe have been classified as SN IIL or SN Iln events, a "unified" model is sought in order 
to explain their properties; at least the properties of their CSM environments and ideally of 
their progenitors as well. 

One of the interesting aspects of the SLSN events that seem to be surrounded by dense, 
optically-thick shells is that the nature of the underlying supernova is hidden and thus re- 
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mains obscure. There seems to be too much mass to be associated with any underlying 
white dwarf thermonuclear explosion, but that leaves many unanswered questions. What 
are the progenitor mass and radius, the ejecta composition and energy? What is the explo- 
sion mechanism, core-collapse or something else? The collective enigmatic nature of these 
super-luminous events means that we have to return to rudimentary studies to explore the 
parameter space that may be appropriate. For this reason we have generalized the basic 
light curve models of A80, 82 using the first law of thermodynamics and the diffusion ap- 
proximation. Where Arnett first proposed his models in the context of radioactive decay, 
we have adapted his technique for an arbitrary prescription of the power input. This allows 
us to explore radioactive decay, shell-shocks, magnetars and, within some limitations, com- 
binations of these power inputs. These basic models allow us to address questions such as: 
what is needed to make the light curve rise have a certain shape, and, independently, what 
shapes the decline; what is needed to generate a nearly symmetric light curve, or a very 
asymmetric one; what is the effect of the initial radius of the supernova progenitor; what are 
the constraints on the density profiles in the supernova ejecta and the circumstellar matter? 
The most general question we propose to explore is whether or not the heterogeneity of the 
current sample of SLSN requires very different progenitors and physics, or if there is some 
common theme expressed in different ways. The goal of this study is not to compete with 
more realistic radiation hydrodynamics simulations but rather complement them by provid- 
ing approximate solutions that will help us understand the importance of the parameters 
involved in the case of SLSNe. Benchmarking against more accurate, numerical results is 
discussed in this work as a means to illustrate the uncertainties and the limitations of our 
model. 

We organize the paper as follows: In §2 we present the general physical assumptions 
of the basic LC model and wc provide solutions in the cases of homologously expanding 
material and fixed SN photosphere. In §3 wc present a variety of physically-motivated 
luminosity inputs and we develop model LCs for shock energy deposition resulting from SN 
ejecta-CSM interaction incorporating the effects of radiative diffusion. We also develop a 
hybrid LC model with ejecta-CSM interaction plus ^^Ni and ^^Co radioactive decay. In §4 
we provide a characteristic fit of the hybrid ejecta-CSM interaction radioactive decay input 
to the LC of the SLSN 2006gy and discuss the implications of this model for the nature of 
the event. Applications to other SLSNe will be presented in a subsequent paper. Finally, 
in §5 we summarize our conclusions. Details of the derivations of the analytical models are 
given in Appendices A and B. 
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2. 



GENERAL LIGHT CURVE MODEL 



Following the prescriptions of A80, 82 for Type I SN LCs, we present a solution for a 
general heating input. The initial assumptions for our model remain the same as in A80, 
82, those being: 1) homologous expansion of the ejecta 2) centrally located power input 
source 3) radiation pressure dominant. Under these assumptions we consider the first law of 
thermodynamics: 



where E = aVT* is the specific internal energy, P = {l/3)aT^V is the pressure, V = p^^ is 
the specific volume where p is density, ej„p is the specific input energy generation rate, L the 
output radiated luminosity and m is the mass coordinate of the fluid element. In general, the 
temperature profile of the diffusion mass, T, ej^p and L are functions of position, x — r/R, 
and time, t, where x is the dimensionless position coordinate for a fluid element relative 
to a flducial radius taken to be the radius of the photosphere. For homologous expansion 
R — Ro + vt where Rq is the initial radius of the ejecta at the moment of shock breakout and 
V is the characteristic expansion velocity of the ejecta. The velocity v is not necessarily the 
photospheric expansion velocity, Vph, as measured from SN spectra, but we can use v — Vph 
as an approximation in some cases. 

For the output luminosity we use the radiation diffusion approximation 



where A = l/np is the mean free path with k. being the mean opacity that we take to 
be a constant, p the density, and c the speed of light. In the following analysis we alter 
the two flrst criteria of ABO, 82 by considering v ^ (flxed photospheric radius instead of 
homologous expansion) and input sources that terminate due to their movement through 
the diffusion mass such as a forward shock that breaks out from an optically-thick CSM 
envelope. The assumption of constant opacity is a weakness of our analytical approach, but 
it was necessary in order to obtain a separable PDE. In more realistic situations that are 
accounted for in radiation hydrodynamics models opacity is depth and time dependent. We 
note that while Moriya et al. (2011) use their radiation hydrodynamics code to compute the 
ionization state and opacities of the underlying explosion, they adopt a Thomson electron 
scattering opacity for fully ionized solar metallicity material {k ~ 0.33 cm^ g^^) within the 
ionization front of the CSM. We will use the same value for the optical opacity throughout 
this work. This general LC model has been considered and expanded by Bhnnikov & Popov 
(1993) for the case of a piecewise constant opacity in power-law density distributions for the 
SN ejecta. As pointed out in their paper, the correct approach to the subject belongs to the 
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class of moving boundary problems and not the usual eigenvalue formulation that has been 
considered so far for the analytical models. A different approach has also been considered by 
Popov (1995) where he adopts a mixed boundary condition for radiation (not the radiative 
zero solution that is considered in the A80, A82, A96 models) and obtains solutions for 
general heating terms in the form of Green's functions but specifically focusing attention on 
radioactive heating. 



2.1. Solution for homologously expanding photosphere 

As shown in Appendix A, the general full solution for the output luminosity from the 
photosphere of the SN ejecta can be written as: 



Lit) = -^e '^i \ e'i < f{t')[— + -]dt' + llf^^e"^* (3) 
td Jo '^td td to 



where we have introduced the effective light curve time scale, td = \/2toth, and where 
to — 3KRo/Vooac — kM/(3cRo is the diffusion time scale, Vqo is the initial central specific 
volume of the ejecta, a is a constant arising from the separation of variables of the LC PDE, 
th — Rq/v is the expansion time-scale, Lq is the initial luminosity input (see Appendix A), 
Eth,o the total SN explosion energy and ^ is a constant that accounts for the density profile 
of the diffusion mass. A80, 82 adopt /3 = 13.8 as a good approximation for a variety of 
diffusion mass density profiles. The term R{t)/Roto (see Equation A9 in Appendix A) can 
be written as (Rq + vt)/Roto or {2Ro/vtd + 2t/td)/td (as we have done in Equation 3) as 
convenient. For Rq ^ {) the additive term Rot/vt^ that appears in Equation 3 vanishes and 
the solution for small initial radius is recovered, as presented in A80, 82. Arnett & Fu (1989) 
and A96 also implicitly present this general solution and apply it to the case of the LC of 
SN 1987A for several power inputs (pulsar heating, swept- up luminosity and fallback). The 
second term in Equation 3 is an initial value term related to the internal energy, Ethfl that 
the SN possesses at a given time. In Equation 3 this term is governed by the energy the 
ejecta have at time t = that they gained from the SN blast wave that subsequently diffuses 
from the optically-thick expanding envelope. In A80, this term is used to explain SN IIP 
LCs, for which the CSM shock input energy from radioactive ^^Ni and ^®Co decay is small 
compared to the total energy of the ejecta at early times. Due to the fast exponential decay 
this "fireball" term may be ignored for our purposes since it affects the LC only at very early 
times. We note that this term does play a role in models for which the input ceases and 
diffusion controls the luminous output. In this context, the coefficient of the exponential 
term becomes Lq, the initial value of the luminosity at the time the power input ceases, after 
which the luminosity decays according to the exponential diffusion term (§3.1, §3.2.1). In §3 
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we will explore several physical luminosity inputs to obtain the final generalized model LCs. 



Next, rather than homologous expansion, we assume v = 0, R = Rq = Rph (and 
thus R{t)/Roto — 1/to) for a fixed photosphere radius. The motivation for this is the fact 
that around massive SN progenitors there can be an optically-thick CSM shell, so that 
when the SN ejecta collide with that surrounding medium, a forward shock is formed that 

propagates into the nearly stationary CSM and a reverse shock propagates into the cjccta, 
both depositing kinetic energy and heating the interacting media. These shocks provide a 
natural source for the output luminosity of the event (Chevalier 1982; Chevalier & Fransson 
1994; Chugai & Danzigcr 1994). Although the diffusion problem that wc solved above 
involves an homologously expanding SN photosphere, we assume the same principles hold 
for a stationary photosphere and diffusion mass. This assumption modifies the final PDE 
that we solve (see Appendix A) and leads to the following general solution: 



where Eth^ is, again, the initial internal energy from the SN blast wave that affects the LC 
at early times and that we ignore in our analysis. For small initial radius {Rq — >■ 0) that 
term also goes to zero. 



Having the general solutions for the output SN LCs provided by Equations 3, and 4 
for homologously expanding and fixed photosphere, respectively, we now consider several 
physical heat inputs as the "source functions" for the integrals of these expressions in order 
to obtain the final output LC models. From the form of these general solutions it follows 
that, in order to obtain a physical LC that rises smoothly to maximum and subsequently 
declines, the luminosity deposition function must be a smooth continuous function that does 
not monotonically increase. For a constant or monotonically rising power input, the output 
luminosity will always increase with time with a monotonically increasing slope until the 
heat input vanishes at some point. For a monotonically declining input, the luminosity will 
increase to a maximum value and then decline, dominated by diffusion and cooling. As shown 
in §3.3, plausible input sources can rise and then decline. At very late times, the output 
luminosity will be the same as the input luminosity as the diffusion time becomes short 



2.2. Solution for fixed photosphere 





3. 



POWER INPUTS 
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compared to the total elapsed time. This property is used as a diagnostic for radioactive 
decay of ^^Ni and ^^Co for some SN LCs where late-time photometric observations indicate 
a decline rate consistent with that of ^^Co (Colgate & McKee 1969; Colgate et al. 1980). 
The same is expected to be the case for other continuous power inputs. 

We may consider two general categories of power inputs: centrally located ones and 
moving ones. A well-known input is the radioactive decay of ^^Ni and ^^Co that is believed 
to power the SN la, SN Ib/c and some SN II LCs (A80, A82), particularly the maximum 
of SN 1987A and SN lib. Although there can be some outward mixing, the newly-formed 
nickel is often taken to be confined near the center of the SN cjccta, around the core of the 
progenitor star. This assumption agrees well with some simulations of SN la, SN Ib/c and 
SN II events (Nomoto et al. 1984). Another centrally located input that has been considered 
for SNe is radiation from a magnetic dipole associated with a pulsar (Ostriker & Gunn 1971) 
as applied in the case of SN 1987 by Arnett & Fu (1989) or a magnetar as proposed recently 
as the power source for the SLSNe 2007bi and 2008es (Kasen & Bildsten 2010; Woosley 2010; 
see also Maeda et al. 2007). 

Moving inputs include the energy released in SN ejecta due to the recombination of H 
and He and shock heating. Arnett & Pu (1989) presented a semi-analytical model that takes 
the effects of H recombination into account in order to reproduce the LC of SN 1987A. The 
recombination front recedes into the expanding SN cjccta. The effect of that is the creation 
of a plateau phase in the SN LC that is more pronounced for large initial radius and large 
mass of the H envelope. The analytical model presented in Arnett & Fu (1989) can therefore 
be used to reproduce some SN IIP LC characteristics. Analytical models for SN IIP LCs 
have also been provided by Popov (1993) and Kasen & Woosley (2009). 

Shock heating can be another important moving power source for some SN LCs. SN 
Iln show evidence for CSM interaction in their spectra with prominent H and sometimes He 
emission features indicating the presence of a shock running into the CSM and depositing 
kinetic energy (Chevalier 1982; Chevaher & Pransson 1994). A significant fraction of the 
recently discovered SLSNe are SN Iln and show signs of CSM interaction, so an analytical 
model that describes the output light curve that results from this power input might be help- 
ful in understanding the basic properties of the CSM involved in the process. A peculiarity 
of this input is that its dynamics depend on the physical characteristics of the interacting 
media (the SN ejecta and the CSM shell). In addition, it is not a continuous input; once the 
reverse shock sweeps up all the available ejecta mass or the forward shock sweeps up all the 
available CSM mass there is no further heating and all that is left is just the diffusion of the 
previously shock-deposited energy. Incorporating forward shock heating into an optically- 
thick CSM envelope requires a different treatment in order to account for the movement of 
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the input source, and the resulting PDE is not formally separable. For the scope of this 
work, we will assume that reverse and forward shock heating are both centrally located, but 
that they terminate; forward shock heating terminates when the shock breaks out of the 
CSM and reverse shock heating when the available SN ejecta mass has been swept-up. With 
these assumptions, the PDE becomes separable (Appendix B). The assumption of centrally 
located power source for the case of the forward and the reverse shocks, although convenient, 
is not generally true and thus increases the uncertainties and limitations of this approximate 
model. For this reason in §3.3.4 we compare our results with numerical results presented by 
Chugai et al. (2004), Woosley, Bhnnikov & Heger (2007) and Moriya et al. (2011) for the 
same initial conditions and discuss the differences. These approximate solutions give useful 
guidance, as we will see, but a proper solution calls for numerical radiation hydrodynamics. 
A combination of power sources, centrally located and moving/terminated, is probably in- 
volved for most core-collapse SN. We consider such a case by developing a hybrid radioactive 
decay and CSM interaction diffusion model in §3.3.3. 



3.1. Instantaneous luminosity input at the SN photosphere 

We start by considering the simple initial value problem where the initial luminosity Lq 
is introduced at t = t„iax at the SN photosphere. This problem was also considered by A80 
and deals with the resulting SN LC if the energy input were just the initial SN shock energy 
for which Lq = Eo/to in Equation 3. This input produces a light curve with instantaneous 
rise to maximum light at t = tmax followed by simple diffusion. A96 also presented solutions 
of this "expanding fireball" problem that has the following form in the homologous expansion 
case: 

r (t—tmax) , 2JIq (t—tmax) i 

L{t) = Loe ^ *d \ (5) 

In the fixed photosphere case, this reduces to: 

L(i) = Loe^^*-*"^"-)/*". (6) 

These results are the same in form as the second terms of Equations 3 and 4 for homologous 
expansion and fixed photospheric radius, respectively. Pure diffusion on the decline has the 
same form for any luminosity input that vanishes after some time interval. Similar diffusive 
decay from an initial value of luminosity at a given time pertains to models in which the 
power input is truncated (§3.2.1, §3.3.2). This simple solution has been considered by some 
authors as an interpretation for the dechne of some SLSNe (for an example, see Smith & 
McCray 2007 for the LC of SN 2006gy). Model LCs powered by initial value luminosity 
input are shown in Figure 1 for the cases of homologously expanding matter and stationary 



photospheres. This model does not provide a natural explanation for a LC with a smooth 
extended rise. 



3.2. Centrally located power inputs 

Centrally located power inputs such as the radioactive decay of ^^Ni and ^^Co as well 
as magnetar spin-down have been considered in previous work. In the following subsections 
we summarize those past results in order to use them for comparison with the models that 
we will introduce. 



3.2.1. Terminated constant CS shock luminosity 

Motivated by the fact that an instantaneous shock luminosity input does not provide an 
output SN LC with a smooth rise to maximum light, we consider a "top-hat" model in which 
a constant input is provided for a finite time and then shuts off. To be specific, we consider 
the input to be a shock with constant luminosity, L^h lasting for time, tg^, and thus producing 
total energy Egh = Lshtsh- Therefore, the input shock luminosity that is provided to the 
expanding diffusion mass is Lsh = Esh/tsh for t < tgh and Lgh = otherwise. This is the 
input luminosity function that wc insert into Equations 3 and 4 to calculate the output LC 
in the cases of homologously expanding matter and stationary matter with a photosphere 
of fixed radius, respectively. We also assume that the shock luminosity deposition takes 
place deep within the diffusion mass, therefore we neglect any corrections attributable to the 
movement of the shock towards the photosphere. Direct integration yields 

in the case of homologously expanding matter and 

Lm-JS^f'- '^'-'^ (s) 

Ht) - <^ |^g-t/to[e*..Ao _ 1], t> tsh, ^ ^ 

\ ^sh 

in the case of fixed photospheric radius. Examples of model LCs for those two cases within 
the context of this type of luminosity input are shown in Figure 2. 

In this type of model the energy that powers the SN hght curve is produced by the 
diffusion of shock-generated energy through an optically-thick CSM shell of large initial 
radius. Smith & McCray (2007) adopted a Loc rise for their model of SN 2006gy based on 
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the early portion of the diffusion models of A82. They did not self-consistently consider the 
input necessary to drive such a rise. This constant power model is related to the "top hat" 
magnet ar- input model that is considered in Kasen & Bildsten (2010) with the exception 
that we solve for the general case of large initial radius and the diffusion time is defined 
somewhat differently (see Appendix A). We note that in the case where Rq is small, the 
Kasen & Bildsten result for small initial radius is recovered. While not especially realistic, 
this "top-hat" model captures the essence of the shell-shock model on both the rise and 
dechne. 



3.2.2. Radioactive decays of^^Ni and ^^Co 

Diffusion of the radioactive decay energy of newly synthesized ^^Ni and ^®Co has long 
been considered the source for powering the LCs of SN la, SN Ib/c and some SN II. A79, 
A80, A82 were the first to analytically solve this problem and later works by Valenti et al. 
(2008) and Chatzopoulos, Wheeler & Vinko (2009) provide the generalized mathematical 
expressions for the output LC including the contribution due to nickel decay and positron 
heating (assuming the same deposition function) as well as the effects of gamma-ray leakage. 
The output luminosity in this case is found to be: 

td Jo vtd td 

j,eco / [— + -]e'^^ e-*'/*c.^t'](i _ ^-M--)^ (g) 
Jo vtd td 

here Rq is the initial radius of the progenitor, Mjsn is the initial nickel mass, tj^i — 8.8 days, 
tco — 111.3 days, e^i — 3.9 x 10^° erg and eco — 6.8 x 10^ erg are the energy 

generation rates due to Ni and Co decays respectively (Valenti et al. 2008). The factor 
(1 — e~^^ ^) accounts for the gamma-ray leakage, where large A means that practically all 
gamma rays and positrons are trapped. The gamma-ray optical depth of the ejecta is taken 
to be r-y = KjpR = At~'^, where is the gamma-ray opacity of the SN ejecta. 

Within this model, for which the time-scales of the power inputs are known, the following 
expression for the mass of the SN ejecta can be obtained: 

where Mej is the mass of the SN ejecta, and /3 is an integration constant equal to about 
13.8. This equation for the diffusion mass (the mass of the SN ejecta in this case) holds for 
both homologously expanding matter and fixed photosphere, but only under the assumption 
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of a centrally located power input. If the power input source moves within the diffusion 
mass, to is itself time-dependent and the diffusion mass also changes with time, depending 
on the position of the source relative to the that of the photosphere. Although this model 
can provide good formal fits to observed LCs of some SLSNe, it cannot be adopted as a 
general explanation due to the fact that in many cases the derived Mjvi is greater than Mej 
for reasonable choices of k. 

The usual assumption of a radioactive decay diffusion model is homologous expansion 
of the SN ejecta. We alter this criterion and consider also the solution for a LC powered 
by the radioactive decays of ^^Ni and ^^Co, but for the case of fixed photospheric radius. 
This may be the case for some SN Iln progenitors that are surrounded by a dense extended 
optically-thick CSM envelope. In such case, the photosphere of the relative diffusion mass 
is coincident with the photosphere of the optically-thick CSM which, in principle, does not 
expand homologously or moves with a velocity small compared to the velocity of the SN 
ejecta. Although those types of events are primarily powered by shock deposited energy 
there may still be some contribution from radioactive decay. In this case, we make use 
of Equation 4 for a radioactive ^^Ni and ^^Co decay input and we arrive at the following 
solution: 

L{t) = ^e'Mie^, - eco) f ^e^e-*'/*-^^' 
Jo h 



3.2.3. Magnetar spin-down 

Recently, Kasen & Bildsten (2010) and Woosley (2010) considered models of SLSN LCs 
powered by the spin-down of a young magnetar (see also Ostriker & Gunn 1971; Arnett & 
Fu 1979; Maeda et al. 2007). In such a model, the energy input by the magnetar is given 
by the dipole spin-down formula: 

E / — 1 

^^"^^^^ ^ til +t/t,r ^^^^ 

where Ep is the initial magnetar rotational energy, tp is the characteristic time scale for spin- 
down that depends on the strength of the magnetic field and / = 2 for a magnetic dipole. For 
a fiducial moment of inertia, the initial period of the magnetar in units of 10 ms is given by 
Pio = (2 X 10™erg/s /Ep)^-^. The magnetic field of the magnetar can be estimated from Pio 
and tp as Bu — (1.3Pio/^p,2/t-)°'^, where Bu is the magnetic field in units of 10^^ G and tp^yr 
is the characteristic time scale for spin-down in units of years. Incorporating the magnetar 
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spin-down deposition function presented here into Equation 3 and including the effects of 
large initial radius, we arrive at the following solution: 

Lit) ^ -e *d "^d I e "'d^ \- z -ro^;, (13) 

where x = t/td and y = td/tp with td again being an "effective" diffusion time. As was the 
case for the radioactive decay diffusion model, the mass of the SN ejecta for this model is 
also given by Equation 10. 

Figure 3 shows an example of a radioactive decay diffusion model LC (solid black curve) 
compared with a magnetar spin-down model LC (solid red curve) where the parameters 
have been chosen so that the output LCs have approximately the same peak luminosity. 
The parameters that were used for the models presented in Figure 3 are Mjvi = 10 M©, 
Ep — 1.37 X 10^^ erg, tp — 30 days and for the homologously expanding case for to — 40 days, 
Rq — 10^^ cm and v — 10,000 km s~^. Whereas the decay times of ^^Ni and ^^Co are known 
experimentally, the magnetar spin-down model contains two adjustable time-scales, tp and 
td, instead of just td as is the case for the radioactive decay diffusion model. The magnetar 
model therefore provides more freedom for fitting to observed SN LCs. 

3.3. Shock heating from CSM-ejecta interaction 

SN Iln and the recently discovered SLSNe show a variety of LC characteristics in terms 
of maximum luminosity, duration, shape and decline rate. Motivated by that, we attempt to 
use the ejecta-CSM interaction scenario as introduced by Chevalier (1982) and Chevalier & 
Fransson (1994) coupled with diffusion as treated by A80, 82 in order to obtain model LCs 
for this type of events. 

3.3.1. Forward and reverse shock luminosity from ejecta-CSM interaction 

We adopt the scenario introduced by Chevalier & Fransson (1994) in which the progen- 
itor star is embedded in a CSM shell described by a power-law density profile pcsM = qr^"^, 
where pcsM is the density of the CSM medium, g is a scaling constant and s the power- 
law exponent for the CSM density profile. In general, q = pc5M,i^i, where pcsM,i is the 
density of the CSM shell at r = ri. We use as a fiducial value for ri = Rp where Rp the 
radius of the progenitor star. Thus we set the density scale of the CSM, pcsM,i, immedi- 
ately outside the stellar envelope. For s = 2, a steady-wind CSM model is recovered, where 
q — M/^AnVyj) with M the pre-SN wind mass- loss rate and Vyj the pre-SN wind velocity. 
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Deviations from s — 2 may indicate a different pre-SN mass-loss history resulting in CSM 
clumps or shells. Values of s close to zero may indicate the presence of CSM "bubbles" or 
shells formed by strong stellar winds as proposed for the CSM environments around Wolf- 
Rayet stars (Chevalier & Liang 1989; Dwarkadas 2011). The density profile of the SN ejecta 
is taken to have double power-law profile where the outer ejecta have a power-law density 
profile, psN = g'^t"'~^r~"' , where g"' is a scaling parameter for the ejecta density profile, 
= l/{4n{6 - n))[2(5 -6){n- 5)Es^]("-3)/7[(3 -6){n- 3)Mej]^''~^^/\ n is the power-law 
exponent of the outer component, and 5 is the slope of the inner density profile of the ejecta 
(values of 5 = 0, 2 are typical), Esn is the total SN energy and Mej is the total SN ejecta 
mass. This assumption for the evolution and the density profile of the SN ejecta is supported 
by realistic numerical calculations. The parameter n varies depending on the nature of the 
progenitor and the presence or absence of a convection zone in the outer parts of the star. A 
fiducial value for n is 11.7 corresponding to the case for red supergiant progenitors (Matzner 
& McKee 1999) whereas lower values of n correspond to the envelopes of more compact 
progenitors. 

The interaction between those two media, the SN ejecta and the CSM, each with power- 
law density profiles results in a forward/circumstellar and a reverse/ejecta shock, the dynam- 
ics of which are described by self-similar solutions presented in Chevalier & Fransson (1994). 
We use those similarity solutions to derive the following expression for the luminosity that is 
produced from this process (neglecting the second, initial value, "fireball" term in Equation 
3; see Appendix B for the full derivation): 

, . ZiTT 5— s n— 5 , V n , x c; „ . 5 — s , ^ 2n+6s— ns — 15 

Linpit) = g'^—sq—sin - 3)^(n - 5)/3|-M^(i + U) e{tFs,* - t) + 

[n s) 

„ , , 5 — n ^K^^/S — 5,0/ V 2n-\-6s — ns — 15 / ^ .\ 

27r(^) — + k) — eitRs,* - t), (14) 

where Pf, and A are constants that depend on the values of n and s and, for a variety of 
values, are given in Table 1 of Chevalier (1982), 9{tFS,* — t), 0{tRs,* — t) denote the Heaviside 
step function that controls the termination of the forward and reverse shock respectively 
(tps,* and tRs^* are the termination time scales for the two shocks) and ti ~ Rp/vsw is the 
initial time of the CSM interaction that sets the initial value for the luminosity produced 
by shocks where vsn = [10(n — b)EsN /^{n — 3)Mej]2/a;o is the characteristic velocity of 
the SN ejecta, where Xq = rQ{t)/RsN{t) is the dimensionless radius of the break in the SN 
ejecta density profile from the inner fiat component (described by 6) to the outer, steeper 
component (described by n) which is at radius ro{t). The first and the second terms in 
Equation 14 refer to the forward and reverse shock luminosity input, respectively. The 
forward shock termination time-scale, tFs,*, is given by the following expression, assuming 
that the input from the forward shock terminates when all the available CSM has been swept 
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up: 

— o)a'^3-n)/n~s)tA nl{s-3)/{n-s) n-s 

tFS, = I , J f ^ I^^^mS^, (15) 

where Mcsm is the total mass of the CSM. Once the forward shock breaks out from the 
optically-thick part of the CSM envelope, the luminosity deposition from it is primarily in 
the UV/X-ray region of the spectrum, not the optical. The time tps,* can represent the 
breakout time if Mcsm is taken to be the mass of the optically-thick CSM, rather than the 
total mass. We adopt this assumption below. Assuming that the reverse shock terminates 
when all the ejecta are swept-up, the reverse shock termination time-scale, tRs,*-, is given by 
the following expression: 

vsN ,^ _ (3 - n)Me 

where Mg^ is the total ejecta mass. 



^RS,* = ^ (1 - \_..3-n ' )^-"]-^- (16) 



To verify the analytical result of Equation 14, we compare it to the LC that is calculated 
by a two-dimensional simulation of the collision of SN ejecta with a wind CSM component 
(s = 2) as presented by Van Marie et al. (2010). We adopt the same parameters for the 
CSM and the SN ejecta components as in their model OOl {n — 11.7, s = 2, Vyj = 200 km s~^ 
and M — 10~^ Mq yr'^) and present the results in Figure 4. The agreement between our 
analytical result (red solid curve) and their simulation (filled circles) is remarkable given that 
they have used a different technique to calculate the radiated luminosity in their simulation. 



3.3.2. Ejecta- CSM interaction with diffusion 

The spectrum produced by optically-thin ejecta-CSM interaction, is expected to be 
a hard spectrum populated with emission lines produced by forward and reverse shocks 
(Chevalier & Pransson 1994; Nymark, Chandra & Pransson 2009). Spectra with these char- 
acteristics have not been seen in SLSNe so far, so it is natural to assume that the effects of 
radiative diffusion of the shock-generated luminosity are important in explaining both the 
currently observed optical light curves and optical spectra of these events. For this reason we 
develop an analytical model that couples the effects of a CSM shock heat input as derived in 
the previous section, with the prescription of radiative diffusion as treated by A80, 82. Wc 
will consider the case of a stationary CSM photosphere for this problem and assume that the 
photosophere is somewhere in the CSM envelope where it will not, in principle, be rapidly 
expanding. In these models we ignore "backwarming" of the photosphere by any shock that 
has propagated beyond the CSM photosphere. We include the effects of terminated shock 
luminosity input. 
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We start by assuming that the luminosity input from the forward and the reverse shock 
happens deep within the fixed photosphere of the CSM and therefore we treat it as centrally 
located. This assumption coupled with the nature of the CSM we explore means that the 
typical shock crossing time scale, Rgh/vsh, is larger than the effective radiation diffusion time 
scale td- This statement is equivalent to the condition that the CSM optical depth is smaller 
than the characteristic optical thickness of a radiation-dominated and radiation-mediated 
shock, c/vsh, and thus in the regime considered by Nakar & Sari (2010) who calculate shock 
breakout LCs for a variety of SN progenitors (see also Ensman & Burrows 1992). Chevaher & 
Irwin (2011) investigated situations in the opposite regime, Rsh/vgh < td, for a steady state 
wind (s = 2), considering the effect of the propagation of a radiation-dominated forward 
shock into the mass-loss region for the cases where the characteristic radiation breakout 
radius {Rd, as defined in their work) is either smaller or larger than the wind termination 
radius (/?«,). In Chevalier & Irwin (2011), diffusion occurs when the forward shock reaches 
Rd and Rsh/vsh ~ td- Subsequently, the diffusion would proceed more rapidly than the 
shock. This is the regime we consider. Wc note that when the condition Rgh/vsh > td is 
satisfied, the shock is no longer radiation mediated. Ofek et al. (2010) used a related model 
to explain the LC of the Iln SN PTF 09uj. Balberg & Loeb (2011) presented a similar model, 
but for a less dense wind in the context of observational signatures of the UV/X-ray shock 
breakout LC. 

Both the forward and the reverse shocks move through the diffusion mass affecting 
the radiation diffusion time and changing the form of the output LC. Accounting for the 
movement of the sources makes the resulting PDE unseparable and the problem hard to solve 
analytically, and we neglect this aspect in the current models. Wc do, however, account for 
the shock propagation by terminating the shock luminosity input after the specific time 
scales that are given by Equations 15 and 16. We assume that the forward shock input 
within the optically thick part of the CSM terminates at a time given by Equation 15 for 
McsM = McsM,th, where McsM,th is the mass of the optically-thick part of the CSM. This 
shock termination time is close to the time of forward shock breakout, tFs,BO, which formally 
occurs when the shock reaches optical depth tau — Vgh/c. The mass of the optically-thick 
CSM is given by: 



to determine the radius of the photosphere, Rph^ where k is the optical opacity and Rcsm 




(17) 



and we use optical depth. 




(18) 
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is the total radius of the CSM determined by the following formula: 

rRcsM 

McsM = 47rg / r'^'^dr, (19) 

where Mqsm is the total mass of the CSM envelope that we use as the basic parameter. 
Solving Equation 19 with respect to Rcsm, Equation 18 for Rph and setting the results 
back into Equation 17 we obtain the final result for McsM,th as an expression that depends 
on the basic model parameters Mqsm, Q, s, Rp and k. Once the forward shock breaks 
out, the bolometric luminosity will still invoke the full forward/reverse shock input given 
by Equation 14, but the optical luminosity will decline since the bulk of the energy will be 
emitted in the UV/X-ray region of the spectrum. A portion of the optically-thin forward 
shock contribution is expected to be in the optical due to re-radiation and electron scattering 
(Chevalier & Pransson 1994), but we neglect this effect. With this model we can also estimate 
the time and intensity of the rise of UV/X-ray flux due to forward shock breakout. 

Now we can implement the SN ejecta-CSM interaction luminosity input from the forward 
and the reverse shock given by Equation 14 into the first term of Equation 4 that accounts 
for the diffusion through an optically-thick CSM with a fixed photosphere and obtain the 
following final expression for the output model LC: 

L(t) = -e~^ / e^[- -g''^q—^(n-3?(n-5)^i-'A^(t' + ti) 

^0 Jo [n - sy 

x0{tFS,BO - t') + 27r(^)t^/3^-V(^)'(t' + t.) '"^"-"" ^(W,. - t')]dt', (20) 

q n — s 

where, in the case of an optically-thick CSM, — KMcsM,th/ ^cRph- After the times of the 
termination of the shock power inputs, the luminosity decays according to Equation 6 where 
the initial value of the luminosity is now that for the corresponding terms in Equation 20 at 
tFS,BO and tns,*, respectively. 

The analytic model described by Equation 20 for the output LCs for SNe powered by 
SN ejecta-CSM interaction depends on the properties of the progenitor star {Esn, M^j, Rp, 
S and n) and the properties of the CSM (s, pcsm,i, k- and Mcsm)- The fact that there 
are many unknown parameters involved in this problem follows from the complex nature 
of CSM interaction due to the large variety of possible CSM environments (winds, shells of 
any density scale) and the large variety of possible SN progenitors (red or blue supergiant 
stars). This natural variety can explain certain differences in the observed optical LC shape, 
duration and luminosity of SN Iln. In the same context, the faster post-maximum decline of 
a SN IIL LC can be the result of particular conditions in this large parameter space involved 
in SN ejecta-CSM interaction, as we discuss next. 
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3.3.3. Ejecta-CSM interaction and radioactive decay input with diffusion 

Here we consider a hybrid model in which the luminosity input is provided by both 
SN ejecta-CSM interaction and by the radioactive decays of ^^Ni and ^^Co. To find the 
combined luminosity input we add the radioactive decay luminosity input to the forward 
and reverse shock luminosity input given by Equation 14. Figure 5 shows an example of 
this hybrid luminosity input for collision of SN ejecta {Esn — 1-5 x 10^^ erg, Mej — 20 M©, 
Rp — 10^^ cm, n — 12) with an optically-thick massive CSM shell (s = 0, Mcsm — 1 Mq) 
in the case where 1 A/q of ^^Ni is formed in the explosion, which might be the case for a 
hypernova or a pair-instability SN (PISN). The early behavior of the input is controlled by 
the relative contribution between the shock inputs and the radioactive decay, and the very 
late-time decay rate is equal to the ^®Co radioactive decay rate. In the cases where CSM 
interaction is dominant, the early decline of the input simply scales as a power law function 
of time with the power being a function of the values of n and s as given by Equation 14. 

To obtain a model bolometric LC for this hybrid input that can be used for fits to 
observed SN LCs we assume that the radioactive decay deposition takes place within the 
whole diffusion mass {M^j + McsM,th)- We assume the CSM interaction luminosity input 
takes place just within McsM,th, since in this context we assume the shocks to be "frozen" at 
the interaction region. We take the interaction region to be the interface between the edge of 
the progenitor star and the CSM envelope and deep within the CSM photosphere. Therefore, 
the final hybrid CSM interaction plus radioactive decay LC model has the following form: 

Lit) = -e *o / e*o [- g-—q—in - ^)\n - b)l3l-'A—{t' + t,) Oitps^BO - t') 

to Jo [i^ ~ 

+27r(^)l^/3rV(^)^(t' + t.f-^'^'^eitns,. - t')]dt' + feim,, 
q n — s '^0 Jo 

xiiem - eco)e-*'/*-' + ecoS-'' /'"^"W , (21) 

where to and t^ correspond to the diffusion time-scales through Mcsm and Mgj + McsM,th 
respectively. Once again, the forward and reverse shock terms are replaced by pure diffusion 
decay terms analogous to Equation 6 when shock input terminates. Figures 6 and 7 show 
some model SN LCs that result from SN ejecta-CSM interaction with a shell (s = 0) and 
a steady-state wind (s = 2), respectively. Examples are plotted for a variety of input 
luminosity combinations so that the effect on the final output LC can be illustrated. The 
parameters used for the models plotted are listed in the captions of the figures. Some of 
the peaks of the resulting model SN LCs are not smooth due to the fact that we have used 
a simplistic Heaviside function prescription for the termination of the forward and reverse 
shock luminosity input. In reality this termination process will be a smoother function of 
time, thus making the peak of those LCs smoother, too. 
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Note in Figures 6 and 7 that when the shock inputs rise monotonically with time (shell 
model; Figure 6) the LC on the rise has a monotonically increasing slope, but when the 
shock inputs decline monotonically with time (wind model; Figure 7), the LC on the rise 
has a monotonically decreasing slope, reminiscent of most observed SNe. The shape of the 
rising LC is thus a potentially strong constraint on the models. Given the large number 
of parameters involved in the hybrid model described by Equation 21, a careful survey of 
the parameter space is needed in order to investigate issues of parameter degeneracy and 
correlation. We will present such an analysis in a follow-up paper that will illustrate fits of 
our model to observed SLSN LCs. From a quahtative perspective, model LCs presented here 
will generally be brighter and briefer for higher values of Esn, Rp and pcsM,i and lower values 
of McsM- Luminosities characteristic of SLSNe are more effectively produced considering 
interactions with dense shells (s = 0) rather than with steady-state winds (s = 2). Models 
with constant density shells tend to be more sensitive to Rp than those with winds because 
of the effect of the forward and reverse shocks. Interactions with steady-state winds are 
capable of producing LCs reminiscent of some normal luminosity SNe Iln, as we will show 
in the follow-up paper. 

We now use the hybrid model presented in Equation 21 to investigate the variety of 
LC shapes that we can obtain from that model and discuss implications for the SN II 
classification scheme. As can be seen in Figure 8, a hybrid CSM interaction plus radioactive 
decay input in which CSM interaction dominates can produce a faster post-maximum decline 
rate that is consistent with the observed rapid post-maximum decline rate ("linear" in the 
logarithmic scale) of SN IIL. This is due to the termination of the dominant luminosity input 
(at t = tns,* in the case illustrated, for which the reverse shock luminosity dominates) that 
will lead to a fast post-maximum decline dominated by cooling and diffusion of previously 
deposited shock luminosity. Furthermore, the differences seen between SN IIL and SN Iln 
in their spectra, namely the absence of P Cygni absorption components in the Balmer lines 
in the latter, maybe due to differences in the optical depth of the CSM involved: A higher 
optical depth would obscure the effects of the underlying SN photosphere expansion thus 
leading to significant CSM interaction and the presence of narrow Balmer emission lines 
without a detectable broad absorption component, whereas a lower optical depth in the CSM 
allows for the underlying SN expansion to be seen and the effects of the CSM interaction are 
mainly imprinted in the LC of the event. Therefore, the SN IIL class may represent a subset 
attributable to a particular restricted range in the large number of parameters involved in 
the more generic SN ejecta-CSM interaction scenario and thus SN IIL may be a subclass of 
SN Iln. 

The generic model presented here might also be used to explain the existence of ap- 
proximately symmetric SN LCs like the ones observed for the transients SCP06F6 (Barbary 
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et al. 2009), PTF09cwl, PTFOQcnd and PTFOQatu (Quimby et al. 2011). An example of 
such a "symmetric" LC around the peak is also seen in Figure 8 when CSM interaction and 
radioactive decay are comparable. This topic will be discussed in more detail in the follow-up 
paper. 

This hybrid model also allows us to predict when the forward shock will break out of 

the CSM photosphere based on the fit to the optical light curve. The shock breakout will 
be followed by an X-ray/UV burst of radiation. The model may also be used to estimate 
the intensity of that radiation and the subsequent UV/X-ray LC. In Figure 9 we illustrate 
this by plotting the optical and UV/X-ray output LCs for a choice of parameters given in 
the captions for the case of SN cjccta interaction with a constant density shell {s = 0; left 
panel) and a wind {s — 2; right panel). 

3.3.4- Comparison with results from radiation hydrodynamics modeling 

The various approximations associated with our analytic hybrid SN LC model and its 
limitations do not allow us to use it to address more specific emission characteristics like radio 
emission, thermalization of the radiation, the ionization state, the change in the opacity of 
the gas, and the evolution of the photosphere properties. Those aspects can be addressed 
only via numerical radiation hydrodynamics simulations and are beyond the scope of this 
project. We can, however, address the limitations and uncertainties of our model bolometric 
SN LCs by benchmarking our results against existing numerical calculations for the same 
initial conditions. To do so, we have used the results from three sets of numerical simulations: 
the Woosley, Blinnikov & Heger (2007) (hereafter WBH07) simulation of the LCs produced 
by a pulsational pair-instability supernova (PPISN), the Moriya et al. (2011) simulations 
of the interaction between SN ejecta from RSG progenitors with optically-thick CSM winds 
and the Chugai et al. (2004) simulation of the collision of SN ejecta with a 0.4 M© CSM 
with density profile pcsM ~ in order to reproduce the LC of SN 1994W. We note that 
in our treatment we neglect the initial offset that the LCs will have in the time axis due to 
the fact that the first fight will emerge after a diffusion time-scale. This offset is shown in 
the simulations of WBH07 and Moriya et al. (2011) and we can also explicitly calculate it 
in each case. 

WBH07 considered the evolution of a 110 Mq star that ends its fife with an oxygen core 
of ~ 50 Mq that is in the domain for which the core is unstable due to electron-positron 
pair production. The collapse and subsequent pulse is not strong enough to disrupt the 
whole star, but only to eject 24.5 Mq of its outer parts. The remainder of the star (with 
radius ~ 10^*^ cm) contracts, encountering pair-instability again after 6.8 years that leads 
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to a second, more energetic pulse (6 x 10^° erg) that ejects 5.1 Mq. The ejecta from this 
second pulse then collides with the slower moving ejecta from the first pulse producing 
a luminous SN-like output LC that was calculated using the code STELLA (Blinnikov & 
Bartunov 1993; Blinnikov et al. 1998; Blinnikov et al. 2006) and presented in supplementary 
Figure 6 of WBH07. The authors consider this PPISN phenomenon as an explanation for 
the SLSN 2006gy. 

To reproduce the WBH07 result, we use our hybrid model, presented by Equation 
21 above, for the same parameters {Esn = 6 x 10^° erg, Rp = lO^'^ cm, M^j = 5.1 Mq, 
McsM = 24.5 Mq, Mn-i = Mq). The parameter pcsM,i is not explicitly given in WBH07 
however we can estimate it, assuming steady-state mass loss. In this case pcsM,i = M/ {iTrVyjRp) 
and M ~ (24.5 Mq)/{6.8 years), so for Vyj ~ 1,000 km s"^ Pcsm,i ^ 1-8 x 10"^^ g cm^^ This 
value is close to the dense shell values seen in supplementary Figures 10 and 11 of WBH07. 
The density profiles of the SN ejecta (second pulse) and the CSM (ejecta from first pulse) 
are also uncertain and not necessarily described by power laws in the context of a PPISN. 
For the purposes of our comparison, we will assume that the SN ejecta from the second pulse 
has a SN-like power law density profile with slope n = 12 and for the CSM we will consider 
two cases: one for a constant density shell (s = 0) and one for a steady state wind (s = 2). 
We also assume in our model that the photosphere is within a stationary CSM while in the 
model of WBH07 the first shell which comprises the CSM is expanding homologously at a 
speed of 100-1,000 km s"^. The results of our comparison are shown in Figure 10. The sohd 
black and red curves correspond to the bolometric LCs for the original ejecta velocity and 
a doubled ejecta velocity, respectively, and the dotted black and red curves correspond to 
the UBVRI LCs as calculated by WBH07. The dashed curves show our analytic LCs for 
s = (green curve), for s = 2 (blue curve) and for s = 2 but for the values of M^j and 
McsM adjusted in order to provide a better fit to the results from the simulations (yellow 
curve). As can be seen, given the uncertainties of the density profiles and the simplifying 
assumptions of our model (centrally located power input and a fixed photosphere in the 
CSM) the two models agree in terms of rise time to maximum light and peak luminosity. 
The post-maximum dechne rates of our models are similar to those found by WBH07, but of 
higher luminosity. A somewhat better agreement is found if we decrease M^j and Mcsm by 
about 50 %, giving some indication of the uncertainty in our ability to estimate parameters. 

Bolometric and color LCs of SNe produced by the interaction of SN ejecta (from a 
RSG progenitor) with a CSM that is a steady state wind (s = 2) or with a s = 1.5 density 
profile slope were presented by Moriya et al. (2011) as results of 1-D numerical radiation 
hydrodynamics calculations. We compare the output LCs from our analytic hybrid model to 
those presented by Moriya et al. (2011) for four of their models: sl3hw2r20m2e3 {Esn = 3 x 
10^1 erg, Rp = bx cm, M^j = 8 Mq, Mcsm = 0.65 Mq, pcsM,i = 2 x lO'^^ g cm'^). 
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Sl5hw2rl0m3e3 {Esn = 3 x 10^^ erg, i?p = 5 x 10^^ cm, Mej = 10 Mq, Mqsm = 0.031 Mq, 
pcsM,i = 2x10-^2 sl5hw2r20m2cl {Esn = 10^^ erg, i?p = SxlO^^ cm, M^j = 10 Mq, 

Mc5M = 0.65 Mq, Pcsm,i = 2 x lO'^^ g cm'^) and sl5hw2r20m3e3 {Esn = 3 x 10^^ erg, 
i?p = 5 X 1013 cm, = 10 Mq, Mcsm = 0.065 Mq, pcsM,i = 2 x lO^^^ g cm'^). The 
details for the parameters used for these models are given in Table 2 of Moriya et al. (2011). 
Figure 11 shows the result of this comparison. The solid curves give the model LCs from 
the Moriya et al. (2011) simulations and the dashed curves our analytical model LCs. Our 
models again produce the rise time to maximum light and peak luminosity as well as the 
width of the main diffusion curve. Our post-maximum decline rates are again consistently 
similar to or slower than the ones found by the simulations. 

Chugai et al. (2004) sucessfully reproduced the plateau-like LC of the SN Iln SN 1994W 
with radiation hydrodynamics simulations that involved contribution from both the radioac- 
tive decay of ^^Ni and interaction with a dense CSM. More specifically they computed three 
models, of which their sn94w58 model provided the best fit to the data. This model con- 
sidered the explosion of a progenitor with a very extended radius {Rp = 1.4 x 10^^ cm; 
their wind termination radius is several times larger) with energy Esn = 1-5 x 10^^ erg and 
Mej = 7 Mq that produced 0.015 Mq of radioactive ^^Ni within a 0.4 Mq CSM with density 
profile slope s = 1. Given those parameters the scale density at the base of the CSM is 
PcsM,i = 1-2 x 10~^^ g cm~^. Using those parameters, we plot our analytical LC model (dot- 
ted curve) and compare it to the Chugai et al. (2004) numerical V-band (solid green curve) 
and U-band (sohd blue curve) LCs in Figure 12. We again see that the overall agreement 
between the models is very good. Our analytical model reproduces the plateau of SN 1994W 
in terms of duration and luminosity. We also observe that the analytical LC model is closer 
to the U-band LC numerical model; it reproduces the two luminosity breaks of the LC as a 
result of forward (early break) and reverse (later break) shock input termination. 

The comparison between our analytic model LCs and radiation hydro LCs for the same 
conditions presented here illustrates the uncertainties involved. The fact that we consider 
the forward and reverse shock power input to be centrally located and the photosphere to be 
fixed in the CSM in order to obtain a separable PDE leads to an overestimate of the diffusion 
time scale in our models because the relative shock and photosphere dynamics in reality will 
lead to an ever decreasing optical depth between the forward shock and the photosphere in 
the CSM. This will have the effect of producing a LC with a slower post-maximum decline 
rate in our models, which is what we see in Figures 10 and 11. Consequently, the overestimate 
of the diffusion time scales means we must employ smaller diffusion masses {Mcsm and Mgj) 
to better fit a given LC. For a given LC, we will tend to underestimate the masses involved. 
This is illustrated by our adjusted s = 2 model to better reproduce the result of WBH07. 
The optical depth and the diffusion characteristics will also be affected by proper treatment 
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of the acceleration of the CSM by the precursor wave, the ionization state of the CSM, and 
the position of the ionization front taking into account the dechne in optical opacity beyond 
it, effects considered by Moriya et al. (2011). 

Our comparison with the models of Moriya et al. (2011) are a particularly stringent 
test, since the CSM in their models was of rather low mass and our models are designed to 
represent the case of especially massive CSM shells. The breaks in luminosity seen in the 
two classes of models are of somewhat different, but closely, related nature. In our case, 
the first break in luminosity occurs when the forward shock breaks out of the optically- 
thick CSM or when the reverse shock has swept up all the available mass of SN ejecta. For 
Moriya et al. (2011), these breaks are the result of the photosphere receding within the 
interaction region; the late time radiation is due to "left-over" thermal emission within the 
ejecta. In our models, the shock moves out to the photosphere; in the models of Moriya et 
al. (2011) the photosphere moves inward to meet the forward shock, thanks in large part to 
the acceleration of the CSM by the radiation flow. For substantially more massive CSM the 
radiative acceleration will be less, and we would expect our models and those of Moriya et 
al. to converge to greater similarity. As noted above, we get quite good agreement with the 
results of Chugai et al. (2004) for models with a more massive CSM. 

Although the physics involved is accounted for in a more accurate way in the numerical 
simulations, our simple analytic results reproduce the basic LC features. We understand the 
sign of the effect of neglecting the decrease in optical depth in front of the forward shock. 
Our models may thus be used as a first step in fitting observed SN LCs to get a basic 
understanding of the parameters involved before proceeding to more expensive numerical 
simulations. 



4. APPLICATION TO THE OPTICAL LC OF SLSN 2006gy 

SLSN 2006gy stirred a lot of discussion among the SN community. It was discovered 
by the ROTSE-IIIb telescope of the Texas Supernova Search project (Smith et al. 2007), 
which obtained unfiltered photometry over the course of ~ 200d providing us with a well- 
constrained LC and explosion date. We allow the explosion date to vary in a limited range 
in our fitting process that will be discussed below. The very late-time decline rate of the LC 
of SLSN 2006gy is reported to be consistent with the decline of ^^Co (Smith et al. 2007). 
A rich database of optical spectra were obtained for SLSN 2006gy (Smith et al. 2007, 2008, 
2010) that provides an extensive record of its spectral evolution. SLSN 2006gy showed strong 
Balmer emission features with their narrow components associated with P Cygni absorption 
indicative of photospheric expansion. The Ho; line profile evolved throughout the course of 
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the LC of SLSN 2006gy showing an evolution that is marked by three phases described in 
Smith et al. (2010). The full width at half-maximum (FWHM) of Ha around maximum 
light reveals characteristic velocities of ~ 4,000 km s~^. We note, however, that this velocity 
information may not directly correspond to the bulk kinematic motion of the SLSN ejecta. 
Smith et al. (2010) discuss the possibilities of line broadening due to electron scattering that 
yields information about the scattering optical depth and thermal motion of the electrons, 
but not about a true physical bulk expansion velocity associated with SLSN 2006gy. We 
will nevertheless use the value of ~ 4,000 km s~^ as a fiducial velocity for our model fitting 
purposes. Spectra of the host galaxy of SLSN 2006gy determined the redshift of the SLSN 
to be z = 0.074, which means that the absolute visual peak magnitude of the event reached 
~ -22 m, making SLSN 2006gy one of the brightest explosions ever discovered. Based on 
the information given by the optical LC and the spectra, SLSN 2006gy was classified as a 
SN Iln event (Smith et al. 2007). 

Various models have been discussed to explain the nature of SLSN 2006gy. Smith & 
McCray (2007) considered a shell-shock diffusion model for which the luminosity output 
around maximum light is reproduced by shock heating due to the interaction between the 
SN ejecta and a dense CSM shell. Their simple model provided a decent fit for the decline 
of the LC, but failed to properly account for the rise of the LC since they considered an 
instantaneous shock input at maximum light. They also considered a contribution from the 
radioactive decay of ^^Co to account for the LC at very late times and estimated a ^^Ni mass 
of about ~ 8 Mq. a pair-instabihty SN (PISN) scenario was discussed for SLSN 2006gy 
in Smith et al. (2007). This model has difficulties in accounting for the large discrepancy 
between the total ejected and nickel mass and it turns out that it is not a good fit to the data, 
as wc will show below. Other models considered interaction with an extensive CSM envelope 
(Smith et al. 2007) or CS clouds (Agnoletto et al. 2009). A popular picture today is that 
SLSN 2006gy resulted from the explosion of a massive LBV-typc star within a dense CSM 
envelope, reminiscent to that of rj Car (Smith et al. 2007). Another model for SLSN 2006gy 
was considered by WBH07 in which the LC of the SN is due to interaction between shells 
ejected as a result of the PPISN process (§3.3.4). We also note that another, exotic model, 
a quark-nova explosion has been considered for SLSN 2006gy and other SLSNe (Ouyed et 
al 2010). 

Given that so far we are lacking a self-consistent LC model for SLSN 2006gy that 
reproduces the whole LC and accounts for its spectral characteristics, we consider the ROTSE 
LC of SLSN 2006gy converted to a pseudo-bolometric LC to test the hybrid model discussed 
in this paper. We assume the bolometric correction BC=0 due to the fact that Smith et 
al. (2007) do not provide such an estimate since there are not adequate simultaneous multi- 
band photometric observations throughout the course of the LC of the event. Therefore we 
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accept the lack of a true observed bolometric LC of SLSN 2006gy as one more uncertainty 
in our fit. We also adopt E{B — V) = 0.72 mag yielding R-band extinction Aji = 1.68 mag 
(Smith et al. 2007). We fit Equation 21 to the pseudo-bolometric LC of SLSN 2006gy and 
present our result in Figure 13. Ideally the fit would be done by a chi-square minimization 
technique (Chatzopoulos, Wheeler & Vinko 2009), but here we have chosen parameters 
by hand to illustrate the capacity of the model to represent the observations. The model 
illustrated was obtained for the following choice of parameters: S — 0, n — 12, s — 0, 
EsN = 4.4 X 10^1 erg, = 40 Mq, Mcsm = 5 {McsM,th = 4.9 Mq), i^^ = 5 x lO^^ cm, 
PcsM,i — 1-5 X lO""-*^^ g cm~^, Mjvi = 2 Mq. For these parameters we estimate the radius of 
the photosphere to be Rph — 2.5 x 10^^ cm and the optical depth of the CSM tcsm ~ 120, 
which is consistent with the fact that the photosphere is within the optically thick shell. 
As argued in the previous section where we presented a comparison with the radiation 
hydrodynamics results of Woosley, Blinnikov & Heger (2007) for the LC of SLSN 2006gy, our 
masses may be underestimated by ~ 50 %. Given the maximum luminosity of SLSN 2006gy 
{Lmax = 1.2 X 10''^ erg s~^) and the radius of the photosphere we derived, we estimate a 
black-body temperature of Tbb — 13,000K. This value for Tbb is consistent with the value 
Tbb — 11,000 - 12,000K that Smith et al. (2007) derived from low-resolution spectra taken 
close to maximum light. We note, however, that black-body emission generally provides a 
poor fit to the SEDs of SNe Iln due to flux dilution and line blanketing and to the fact that 
other emission mechanisms are dominant (see for example Dessart et al. 2009 on SN 1994W, 
Miller et al. 2010 on SN 2008iy; Chatzopoulos et al. 2011 on SN 2008am). 

The fit to the LC of SLSN 2006gy given in Figure 13 was based on a constant density 
shell. To illustrate the utility of our models to limit parameters, we note that it was diffcult 
to provide an equally satisfactory fit based on a finite mass, "steady-state" wind, with s = 2. 
The CSM luminosity input is generally too small to account for the observed luminosities 
even for large values of the CSM scale density, pcsM,i- Extraordinary values for pcsM,i 
(10~^, 10^^, 10"'^ g cm~^) imply a wind that starts at the stellar photosphere, but with a 
higher density than the outer envelope of the star. This condition is unphysical and implies 
extraordinary mass-loss rates (>100 solar masses per year). Even with extraordinary values 
for pcsM,i there are fitting issues because the input for the forward and reverse shocks in the 
s = 2 case declines rapidly monotonically, thus forcing the rise to maximum to be very rapid. 
Figures 7 and 9 (second panel) illustrate s = 2 cases for which the rise time is less than 15-20 
days. For s = 2, it is very difficult to produce a longer rise time without, for instance, having 
a wind mass considerably in excess of the ejecta mass and an overall poor fit in detail. Wind 
models that roughly give the correct maximum luminosity also give a continuously-declining 
CSM shock input that produces a very slow decline of the output luminosity until shock 
termination occurs that does not agree with the observations. We thus find that while the 
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models do have a large number of parameters, the observations considerably constrain the 
model parameters in practice. 

We also plot in Figure 13 (dashed curve) the expected UV/X-ray LC due to the forward 
shock propagation in the optically-thin part of the CSM based on the parameters determined 
by the fit to the optical LC. We see that the forward-shock luminosity, and associated UV/X- 
ray burst, is terminated at about 54 days after explosion in the rest-frame. This may be 
consistent with the low X-ray flux detected by Chandra at about 54 rest-frame days after 
explosion as reported by Smith et al. (2007). Once the forward shock terminates upon 
reaching the outer edge of the shell, the UV/X-ray luminosity falls dramatically. This result 
also implies that the UV/X-ray burst prior to termination of the forward shock would have 
been bright enough to be detectable. 

We conclude that interaction of typical SN ejecta with a massive optically-thick CSM 
shell plus a contribution from 2 Mq of radioactive Ni seems to reproduce the observed LC 
of SLSN 2006gy. In this context, the progenitor of SLSN 2006gy was probably a massive 
star that underwent significant and episodic mass loss. The massive CSM shell could have 
been the result of either LBV-type ejection, or a bubble due to interaction of stellar winds of 
previous epochs or the result of a PPISN. We also note that our hybrid CSM interaction plus 
radioactive decay model fit for SLSN 2006gy is similar to the idea presented by Agnolctto et 
al. (2009) in the sense that they also used arguments for combined CSM interaction (with 
CS clumps in their decription) and radioactive decay of less nickel than previously implied 
(3 Mo). 

5. DISCUSSION AND CONCLUSIONS 

We derived semi-analytical models for the LCs of SN II by incorporating the effects of 
diffusion as treated by A80, 82 with a variety of inputs that arc considered candidate SN LC 
powering mechanisms (radioactive decays of ^^Ni and ^^Co, magnetar spin-down and cjccta- 
CSM interaction). We found solutions for the cases of a fixed (in radius) CSM photosphere 
and homologously expanding matter, including the effects of terminated power input in the 
case of CS shock heating. 

One implication from our results is that the principal factor that gives rise to the 
observed diversity in the shape, duration and luminosity of SN Iln is the presence or absence 
of a CSM enviroment and the physical properties that are associated with it; mainly its 
optical depth, characteristic density and density profile. SN IIP, SN lb and SN Ic are all the 
result of core-collapse explosions with the difference that in SN lb and SN Ic events the H 
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and He envelope of the progenitor star is lost and in SN IIP the large outer H envelope is 
retained. The CSM enviromcnt around these events is probably of very low optical depth, 
so that the effects of ejecta-CSM interaction are negligible and the only mechanism that 
gives rise to the optical luminosity output are shock heating of the cjccta (for SN IIP) plus 
an internal source (radioactive decays or magnetar spin-down) combined with the effects of 
H and He recombination. For SN IIP events where the H envelope is retained, the effects 
of H recombination are more pronounced, while for SN lb and SN Ic events the LCs lack 
an extended plateau phase due to the absence of the outer envelope. Some SN Ib/c events 
develop SN Iln characteristics in their nebular spectra at late times, indicating that the 
ejecta have finally reached the previously expelled H or He envelope and started to interact 
with it (Pastorello et al. 2008). 

On the other hand, in SN IIL and SN Iln events, the CSM envelopes have higher optical 
depth so that the effects of ejecta-CSM interaction dominate the output luminosity. In the 

case of SN IIL the optical depth is probably moderate since P Cygni features are seen and the 
LC decline is faster, while for SN Iln a high optical depth is implied that keeps the expansion 
of the SN ejecta obscured and produces a LC of longer duration. The pronounced effects of 
the optically-thick CSM environment on the LC of luminous SN Iln manifest themselves by 
showing a variety in the LC shapes and durations of these events. This leads us to argue 
that the SN IIL classification could be a byproduct of the diversity of CSM properties and 
not a separate class of SN explosions. The faster decline and apparently symmetric shape of 
the LCs of some SN IIL and SLSNe (SCP06F6; Barbary et al. 2009, PTFOQcwl, PTFOQcnd 
and PTFOQatu; Quimby et al. 2011) may be well reproduced by models of shock heating, 
where the forward and reverse shocks terminate their contribution once they have swept up 
most of the available SN ejecta and CSM mass. If that is true, then strong pre-SN mass loss 
might constitute a significant source of diversity among SNe Iln, producing SLSNe at the 
high mass end, and putting such interacting SNe into a different parameter space compared 
to SNe that do not strongly interact. This point of view was advocated by Blinnikov & 
Bartunov (1993) where they presented radiation-hydrodynamics models of SN IIL in order 
to fit the LC of SN 1979C. 

We applied our hybrid CSM interaction plus radioactive decay model in the case of the 
LC of SLSN 2006gy. We found that the LC of this extraordinary event could be reasonably 
well reproduced by a model where interaction of ~ 40 Mq of SN ejecta with a circumstellar 
shell of ~ 5 Mq plus radioactive decay of ~ 2 M© of ^^Ni provide the luminosity input. 
Although we derive a smaller ^^Ni mass than previous authors, it remains the case that 
SLSN 2006gy was a brilliant explosion and that the progenitor star must have undergone 
episodic mass- loss, a picture consistent with an LBV- type star (Smith et al. 2007). We note 
that alternative scenarios for the nature of the massive CSM shell around SLSN 2006gy, such 
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as a wind blown- bubble (Dwarkadas 2011) or a PPISN (Woosley, Heger & Blinnikov; 2007) 
cannot be excluded. 

The models we present here are only approximate and have a large number of parame- 
ters, but they also provide the means to efficiently explore a large range of parameter space. 
In this way, parameters can be chosen to guide more elaborate and expensive numerical 
radiation-hydrodynamics calculations. In a subsequent paper we will present fits to the LCs 
of many of the recently discovered SLSNe and other peculiar transient events. The ultimate 
goal is to deduce the physical properties of the progenitors and the CSM environments of 
a variety of SNe using their observed optical LCs and spectra in order to constrain the 
mass-loss history and pre-SN evolution of massive stars. 
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A. Derivation of the general LC model for centrally located power sources. 

A80, A96 presented a solution for light curves powered by radioactive decay of ^^Ni 
based on the energy equation: 

dL 

E + PV^ej,i-—, (Al) 

where E is the specific internal energy, P is the pressure, V = p^^ is the specific volume, eNi 
is the specific energy generation rate corresponding to the radioactive decay of ^^Ni and the 
luminosity is given by the diffusion approximation: 

4nr^XcadT'^(x,t) , , 

^ = — 3 h^' '^2' 

where k is an appropriate opacity. A80, 82 solved this set of equations by assuming 1) 
homologous expansion, 2) a centrally-located power source, and 3) that radiation pressure 
was dominant. A80, 82 did not consider a general power input source as we need to consider 
here. 

Kascn & Bildsten (2010) presented a related solution for the hght curve based on an 
integrated version of the energy equation and a general power input, of which a magnetar 
was a specific example. Kasen & Bildsten did not show explicitly how this global, integrated 
solution is related to the solution of A80, 82 that was expressed in terms of local, specific 
quantities. Here we present the most general solution in the formulation of A80, 82, show 
how that solution is generalized to any power input, and how that solution relates to the 
global, integrated solution. 



This preprint was prepared with the A AS lAT^jX macros v5.2. 
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Following A80, 82 we adopt the same prescription for separation of variables, 

T\x, t) = T^oHxm) [i?o/(i?o + vt)]\ (A3) 

where x — r/R is the dimensionless radial variable, Too is the initial central temperature, v 
is the scaling velocity of the homologous expansion v{x) = xv, and Rq is the initial radius. 
Adopting K, = constant, we can then write, following the analysis of ABO, 82, the luminosity 
at the surface, function of time as: 

L(l,t) = f e^i^'^^f{t')[^]dt' + Elh^e-^''/'l-^'^'/<\ (A4) 

Im Jo toRo Meoblth 

where M is the cjccta mass, cq is the amplitude of the specific power input in erg s~^, 
f{t) is the timc-dcpcndence of the power input, Ethfl = 47r/?QaTQQ0(O)/t/i is the initial total 
thermal energy, is the LC timc-scalc which is the geometric mean of the diffusion time- 
scale to = 3/ti?oPoo/ca and the expansion time-scale = Rq/v, therefore = 2toth- The 
factor R{t)/Roto can be written as (i?o +vt)/Roto or {l/td){2Ro/vtd + ^t/td), as convenient 
to evaluate the hmits Rq ^ and v — > 0. The density is assumed to scale as p{x,t) = 
Poor]{x)[Ro/R{t)]^, where poo is the central density at time and the function a appearing 
in the definition of to, 

a— ^ d . x'^ dijj. 
x'^ijj{x) dx 7]{x) dx ' 

is a constant by separation of variables. The parameter h = C,{x)'r]{x)/ilj{x) is assumed to be 
constant where ^(x) represents the radial distribution of the power input, Ith = Jq ip{x)x'^dx, 
and Im = rj{x)x'^dx. A80, 82 show that to = kM/(5cRq where j3 = AiialM/^ — 13.8 for 
a variety of density distributions and that hMIth/Ihi = M^-, the initial amount of nickel 
injected in the radioactive decay model. 

With these definitions and relations we can write for a general power input 

€inp = eo^{x)f{t), (A6) 

and, expressed in terms of luminosity 

-M 



where 



nlVl 

Linp^ ^inpdm^ Lof{t), (A7) 
Jo 



Lo = Meob^, (A8) 

-'M 



or, formally, 

rt2 2^0* 2Rot' 



Jo toRo to 
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This is the form of the output luminosity we present in Equation 3 of the main text. 
Returning to Equation Al, we can integrate the energy equation over mass to find 

AT:R{tfaT^^[^Y4>Ith = Ur^p - L (AlO) 

where L — (dL/dm)dm. Kasen & Bildsten (2010) began with an approximate integrated 
form of the energy equation to write 

-[-7ri?(t) aT ] + -aT AnR = L,„, - ^- (All) 

where the second term on the right hand side is —L. Using the same separation of variables 
for T as in Equation A3, this equation can be written as 

47ri?(t)=^aTot,[-|^]V^ = U^p - L. (A12) 

This equation agrees with Equation AlO only if ^l){x) — SIth, which is only true if ^0 = con- 
stant. This means that the temperature as a function of radius has to be constant. The 
solution of Kasen & Bildsten is thus not strictly speaking self-consistent, since a temperature 
gradient is required and employed to estimate the luminosity. Their solution is neverthe- 
less a useful approach, given the uncertainties involved in modeling light curves driven by 
uncertain phenomena. 



B. Derivation of the SN ejecta-CSM interaction self-similar luminosity input. 

To derive an expression for the luminosities of the forward and reverse shocks we assume 
that all of their kinetic energy converts efficiently to radiation (radiative shock approxima- 
tion) : 

^^^^Jt {l^''"'"^^'^) ^ MsMVsh + (Bl) 

where Mgw is the swept-up mass behind the shock and Vgh the shock velocity. 

In order to estimate a final expression for the luminosity presented in Equation Bl, 
we need to know the dynamics of the forward and reverse shock (radius, velocity and ac- 
celeration) as a function of time. Chevalier (1982) considered the interaction between two 
media with power- law density profiles: the SN ejecta density profile psN — g'^t'^'^r"'^ , where 
g'^ is a scahng parameter for the ejecta density profile, g'^ — l/(47r(5 — n))[2(5 — 5){n — 
5)£;5Ar]^"~^^''V[(3 -S){n- 3)MeJ]("-^)/^ n is the power-law exponent and S is the slope of 
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the inner density profile of the ejecta (values of 5 = 0, 2 are typical) and the CSM density 
profile pcsM — QT^~^, where pcsM is the density of the CSM medium, g is a scaling constant 
and s the power-law exponent. In general q = pcsM,i'i^i where pcsM,i is the density of the 
CSM shell at r = ri. We use as fiducial value ri = Rp, where Rp the radius of the progenitor 
star. Thus we set the density scale of the CSM, pcsM,i, immediately outside the stellar 
envelope. 

Using momentum conservation, Chevalier (1982) and Chevalier & Fransson (2001) found 
the following self-similar solutions for the radii of the forward and the reverse shocks, re- 
spectively, as a function of time: 



and 



Rpit) = Rp + Pf 



-Rii(^) — Rp + Pr 



(B2) 



(B3) 



where (3^, and A are constants that depend on the values of n and s and, for a variety 
of values, are given in Table 1 of Chevaher (1982). More specifically the parameters /3f and 
refer to the ratio of the shock radius to the radius of the contact discontinuity that forms 
as a result of SN cjccta-CSM interaction (/3^ = Ri/R^ and /3r = R2/RC for the forward and 
the reverse shock, respectively, where Rc is the radius of the contact discontinuity given by 
Equation 3 of Chevalier (1982) and Ri and R2 are the radii of the forward and the reverse 
shocks respectively). Using vf,r = dRF,R/dt, vf,r = dPRF,R/dt^ yields the velocity and the 
acceleration of the shock as a function of time. 

Given the shock dynamics derived above, we can calculate the swept-up mass behind 
the forward shock assuming that Rp is smaller than Rpit) for t > 0: 



Rpit) 



47r/3 



3-s 



3-s (n-3)(3-s) 



/ \ 2 1 ^*^yF n—^ rji— — ' ' 

PcsM[r)r dr = — —q^-olAg J"-t " 



Rp 



3-s 



(B4) 



and the swept-up mass behind the reverse shock: 
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where Vsn is the characteristic SN expansion velocity and Rsn = VsNt is the radius of the 
SN photosphere which is assumed to expand homologously. We have again assumed that Rp 
is small. 
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Substituting Equations B2 through B5 and their derivatives in Equation Bl yields the 
final results for the luminosity input from the forward and the reverse shocks (-^f(^) and 
Lji{t) respectively): 

, . ^TT ^ 5 — s n — 5 , , o , , c „ 5— s , ^ 2n+6s— ns — 15 , , , 

Lrit) = _ g^—^q—^{n - 2,f{n - 5)/3|-M^(t + U) e{tFs,BO - t), (B6) 

[n s) 

and 

/ ^ . 5— n „ , O — 5 , Q - , 2n+6s— ns— 15 , . , . 

Ln{t) = 27r(^)^^"(^)^(i + i^^^^^fe,* - t), (B7) 

where Oitps,* — i) and 9{tRs,* — t) denote the heaviside step function that controls the 
termination of the forward and reverse shock, respectively, and ti ~ Rp/vsN is the initial 
time of the CSM interaction which sets the initial value for the luminosity produced. The 
velocity, vsn, in Equation B5 is given by 

[10(n-5)i^W3(n-3)Me,]i 
Vsn = , (BSj 

where xq = ro{t)/ RsN{t) is the dimensionless radius of the break in the SN ejecta density 
profile from the inner fiat component (controlled by 6) to the outer, steeper component 
(controlled by n) which is at radius ro(t). Lji{t) in Equation B7 is independent of vsn and 
hence of Xq. Summing Equations B6 and B7, the total luminosity input from self-similar SN 
ejecta-CSM interaction is 

Linp(t)^LF(t) + LR(t). (B9) 

The forward shock termination time-scale, tps,*, is given by the following expression, assum- 
ing that the input from the forward shock terminates when all the available CSM has been 
swept up: 

_ o^/j(3-n)/n-s)r /I nl(s-3)/(n-s) n-s 

tFS, = 1^^-^^^^ j^^^ l^^^MS^, (BIO) 

where Mcsm the total mass of the CSM. Following the same assumption, the reverse shock 
termination time-scale Irs * is given by the following expression: 



Vsn /, (3 — n)M( 



1 , n- 



After termination, the luminosity of each component decays in a manner analogous to Equa- 
tion 6 in the main text such that 

Lpit) = LF(tF5,so)e-^*-*^«'^°)/*°, (B12) 

and 

Lnit) = L«(t^5,*)e-(*-*««-)/*°. (B13) 
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Fig. 1. — Examples of LC models for instantaneous luminosity input in the cases of homolo- 
gously expanding (solid curve) and fixed (dashed curve) photosphere (see §3.1). The decline 
is determined purely by diffusion. The initial luminosity for this model is Lq = 10^"^ erg s~^ 
introduced at t = 10 days for the choice of diffusion time = 100 days. For the homolo- 
gously expanding photosphere case the initial radius is Rq = 10^^ cm and the photospheric 
velocity v = 10,000 km s~^. For the fixed photosphere, Ro = 10^^ cm=constant. 
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Fig. 2. — Example of LC models for a terminated constant luminosity input ("top-hat" - 
dotted curve) in the cases of homologously expanding matter (solid curve) and fixed photo- 
sphere (dashed curve). See §3.2.1 for details. The models shown are for total input energy 
2 X 10^^ erg that terminates at tgh — 50 days and for the choice of diffusion time to — 100 days. 
For the homologously expanding case the initial radius is Rq — 10^^ cm and the photospheric 
velocity v — 10,000 km s~^. For the fixed photosphere, Rq — 10^^ cm=constant. For this 
choice of parameters the diffusion mass is determined to be 11.6 Mq. 
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Fig. 3. — Example of LC models for ^^Ni and ^^Co radioactive decay input (dashed black 
curve) and for magnetar spin-down input (dashed red curve; the output LC models are 
presented with solid curves of the same corresponding colors) . The radioactive decay diffusion 
model is drawn for M^rj = 10 Mq and the magnetar spin-down model for Ep = 1.37 x 10^^ erg 
and tp = 30 days. The choices of the model parameters were made such that the two model 
LCs have approximately the same maximum luminosity. Both models are drawn for the 
homologously expanding photosphere case for diffusion time to = 40 days, Rq = 10^'^ cm and 
V = 10,000 km s~^. For details on the models see §3.2. 
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Fig. 4. — Comparison of the model LC from the Van Marie et al. (2010) two-dimensional 
optically-thin radiation hydrodynamics simulation of the collision of SN ejecta from a red 
supergiant progenitor with a steady-state wind (s = 2) CSM component (filled circles) 
with the analytic self-similar result for the FS luminosity for the same CSM and SN ejecta 
parameters. The break at about 10 days in the self-similar model represents the time of the 
forward shock termination. 
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Fig. 5. — Example of a hybrid ejecta-CSM interaction plus ^^Ni and ^^Co radioactive decay 
input (solid black curve). The dashed red, green and blue curves represent the forward, 
reverse and radioactive decay luminosity inputs, respectively. The vertical dotted lines 
indicate the termination time for the forward and reverse shock luminosity inputs. This 
hybrid luminosity input model is drawn for the following choice of parameters: 6 = 0, 
n = 12, s = 0, EsN = 1 X 10^1 erg, Mej = 15 Mq, Mcsm = 1 M©, Rp = 1 x 10^^ cm, 
PcsM,i = 5 X 10-13 g cm-3. Mm = 0.05 M©. 
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Fig. 6. — Comparison of SN ejecta-CSM interaction and hybrid SN ejecta-CSM interaction 
plus ^^Ni and ^^Co radioactive decay LC models for the case of collision with a constant 
density CSM shell (n = 12, s = 0). The hybrid forward, reverse and ^®Ni and ^^Co radioactive 
decay model is shown in black and the pure forward and reverse shock ejecta-CSM interaction 
model in red. In each case the dashed curve represents the input luminosity and the solid 
curve the output LC model. The dotted vertical lines indicate the forward shock breakout 
time (essentially equal to the time of termination of forward shock input in the optically-thick 
part of the CSM) and the reverse shock termination time scale. The forward termination 
time scale is very close to its breakout time scale in the case of a CSM constant density 
shell. The models were drawn for the following choice of parameters: 6 = 0, Esn = 10^^ erg, 
Mej = 12 Mo, McsM = 1 Mq, Rp = 10^^ cm, pcsM,i = IQ-^^ g ^^^-3^ = 0.08 Mq. 
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Fig. 7. — Comparison of SN ejecta-CSM interaction and hybrid SN ejecta-CSM interaction 
plus ^^Ni and ^^Co radioactive decay LC models for the case of collision with a steady state 
wind (n = 12, s = 2). The wind is terminated at a finite mass. The hybrid forward, reverse 
and ^^Ni and ^^Co radioactive decay model is shown in black and the pure forward and 
reverse shock ejecta-CSM interaction model in red. In each case the dashed curve represents 
the input luminosity and the solid curve the output LC model. The dotted vertical lines 
indicate the forward shock breakout time (essentially equal to the time of termination of 
forward shock input in the optically-thick part of the CSM) and the reverse shock termination 
time scale. The reverse shock termination time scale is much later, outside of the range 
of this graph. The models were drawn for the following choice of parameters: 6 = 0, 
EsN = 10^1 erg, M^, = 12 M^, Mcsm = 0.1 M©, = 2x10^' cm, pcsM,i = 5x IQ-^^ g cm'^. 
Mm = 0.05 Mq. 
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Fig. 8. — Comparison of hybrid CSM interaction plus radioactive decay LC models for 
a case where CSM interaction in a constant density shell is the dominant contributor to 
the output luminosity (solid curve), a case where CSM interaction and radioactive decay 
have comparable luminosity inputs (dashed curve) and a case where radioactive decay is the 
dominant input luminosity source (dotted curve). The post- maximum decline rate can be 
more rapid for the hybrid models for which CSM interaction dominates or is comparable 
to the radioactive decay input. Such models are consistent with the decay seen in SN IIL. 
For appropriate choices of parameters, the hybrid model can also produce an approximately 
symmetric LC shape around peak luminosity, similar to the observed optical LC of the 
peculiar transient SCP 06F6 (Barbary et al. 2008) as illustrated here by the model for 
which the CSM dominates. For the case where CSM interaction is dominant the following 
parameters were used: (5 = 0, n = 12, s = 0, Esn = 10^^ erg, Mf,j — 10 Mq, Mqsm = 5 Mq, 
Rp — llO-*^^ cm, pcsM,i — 10~^^ g cm~^, M^i — 0.05 Mq. For the case where radioactive 
decay input is comparable to the CSM input the following parameters were used: 5 — 0, 
n = 12, s = 0, Esn = 4 x 10^^ erg, Mej = 40 Mq, Mcsm = 7 Mq, Rp = 10^^ cm, 
PcsM,i = 10~^^ g cm~^, Mjvi = 2 Mq. Finally, for the case where the radioactive decay input 
is dominant the following parameters were used: S = 0, n = 12, s = 0, Esn = 7 x 10^" erg, 
Mej = 21 Mq, Mcsm = 2 Mq, Rp = lO^^ cm, pcsM,i = 6 x lO'^^ g ^^^-3^ ^ q g 
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Fig. 9. — Comparison between the optical (solid curve) and the UV/X-ray (dashed curve) LC 
for a hybrid CSM interaction plus ^^Ni and ^®Co radioactive decay model in the case of SN 
ejecta interaction with a shell (s = 0; left panel) and a wind (s = 2; right panel). The dotted 
vertical lines indicate the times of forward shock breakout (essentially equal to the time of 
termination of forward shock input in the optically-thick part of the CSM), forward shock 
termination (the end of optically-thin X-ray and UV input by the forward shock) and reverse 
shock termination. The s = model corresponds to the following choice of parameters: 
5 = 0, n = 12, EsN = 1-2 X 10^^ erg, M^j = 20 Mq, Mcsm = 2 Mq, i?p = 2 x 10^^ cm, 
PcsM,i = 5 X 10~^^ g cm~^, M^i = 0.05 Mq. The s = 2 model corresponds to the following 
choice of parameters: 5 = 0, n = 12 Esn = 1-2 x 10^^ erg, Mej = 18 Mq, Mcsm = 1 Mq, 
Rp = 10^4 cm, pcSM,l = 2 X lO'^^ g ^^^-3^ _ q 
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Fig. 10. — Comparison of analytical SN LCs with numerical results from Woosley, Blinnikov 
& Heger (2007) for the same initial conditions. The solid black and red curves correspond to 
the bolometric LCs for the original ejecta velocity and a doubled ejecta velocity, respectively, 
and the dotted black and red curves correspond to the UBVRI LCs as calculated by WBH07. 
The dashed curves show our analytical LCs for s = (green curve), for s = 2 (blue curve) 
and for s = 2 but for the values of M^j and Mcsm adjusted in order to provide a better fit 
to the results from the simulations (yellow curve; see text). 
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Fig. 11. — Comparison of analytical SN LCs (dashed curves) with numerical results from 
Moriya et al. (2011) (solid curves) for the same set of initial conditions for four mod- 
els: sl3hw2r20m2e3 (upper left panel), sl5hw2rl0m3e3 (upper right panel), sl5hw2r20m2el 
(lower left panel) and sl5hw2r20m3e3 (lower right panel). See text for model parameters. 
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Fig. 12. — Comparison of analytical SN LC (dashed curve) with the numerical V-band (solid 
green curve) and U-band (solid blue curve) LCs for model sn94w58 of Chugai et al. (2004) 
for SN 1994W. See text for model parameters. 
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Fig. 13. — Fit of a hybrid ejecta-CSM interaction plus ^^Ni and ^^Co radioactive decay LC 
model (solid red curve) to the ROTSE LC of SLSN 2006gy. The parameters for this model 
are the following: 5 = 0, = 12, s = 0, Esn = 4:A x 10^^ erg, Mej = 40 M©, Mcsm = 5 M©, 
Rp = 5 X 10^^ cm, pcsM,i = 1-5 x 10~^^ g cm"^. Mm = 2 Mq. The dashed curve shows 
the expected UV/X-ray LC for this event corresponding to those fitting parameters. The 
dotted vertical lines indicate the times of forward shock breakout (essentially equal to the 
time of termination of forward shock input in the optically-thick part of the CSM), forward 
shock termination (the end of optically-thin X-ray and UV input by the forward shock) and 
reverse shock termination. 



